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In a combined experimental and theoretical effort, we demonstrate a novel type of dipolar sys¬ 
tem made of ultracold bosonic dipolar molecules with large magnetic dipole moments. Our dipolar 
molecules are formed in weakly bound Feshbach molecular states from a sample of strongly mag¬ 
netic bosonic erbium atoms. We show that the ultracold magnetic molecules can carry very large 
dipole moments and we demonstrate how to create and characterize them, and how to change their 
orientation. Finally, we confirm that the relaxation rates of molecules in a quasi-two dimensional 
geometry can be reduced by using the anisotropy of the dipole-dipole interaction and that this 
reduction follows a universal dipolar behavior. 

PACS numbers: 34.50.-s, 31.10.-l-z, 34.50.Cx 


Ultracold dipolar particles are at the heart of very in¬ 
tense research activities, which aim to study the effect 
of interactions that are anisotropic and long range [1, 2]. 
Dipolar quantum phenomena require ultracold gases and 
a strong dipole-dipole interaction (DDI). So far, strongly 
dipolar gases have been obtained using either atoms with 
a large magnetic dipole moment or ground-state polar 
molecules with an electric dipole moment [2]. With both 
systems, many fascinating many-body quantum effects 
have been observed and studied, such as the d-wave col¬ 
lapse of a dipolar Bose-Einstein condensate [3, 4], the de¬ 
formation of the Fermi sphere [5] , and the spin-exchange 
phenomena [6, 7]. 

Here, we introduce a novel kind of strongly dipolar 
particles. These are weakly bound dipolar molecules pro¬ 
duced from a pair of atoms with large magnetic dipole 
moments, such as erbium (Er). The central idea is that 
these molecules can possess a very large magnetic mo¬ 
ment /r up to twice that of atoms (e. g. 14 Bohr magne¬ 
ton, fj,B, for Er 2 ) and have twice the mass of the atoms. 
As a consequence, the degree of “dipolarity” of the mag¬ 
netic molecules is much larger than the one of atoms. 
This can be quantified in terms of the dipolar length 
Qd = /{Airh?) [1], which solely depends on the 

molecular mass m and on /r; h is the Planck constant 
divided by 27r. To give an example, Er 2 with /r = 14 p,B 
has a ad of about 1600 oq, which largely exceeds the typi¬ 
cal values of the s-wave scattering length. Here, oq is the 
Bohr radius. Moreover, in contrast to ground-state het- 
eronuclear molecules, the dipole moment of the magnetic 
molecules does not vanish at zero external (magnetic) 
field, opening the intriguing possibility of investigating 
the physics of unpolarized dipoles. 

In a joined experimental and theoretical effort, we 
study the key aspects of ultracold dipolar Er 2 molecules. 


including the association process, the molecular energy 
spectrum, the magnetic dipole moments, and the scat¬ 
tering properties in both three- (3D) and quasi two- 
dimensional (Q2D) geometries. 

Erbium belongs to the class of strongly magnetic lan¬ 
thanides, which are currently attracting great attention 
in the field of ultracold quantum gases [4, 8-10]. In¬ 
deed, these species exhibit unique interactions. Beside 
the long-range magnetic DDI, these species have both 
an isotropic and an anisotropic contribution in the short 
range van der Waals (vdW) potential. The latter re¬ 
sults from the large non-zero orbital momentum quan¬ 
tum number of the atoms [11, 12]. This manifold leads 
to an extraordinary rich molecular spectrum, reflecting 
itself in a likewise dense spectra of Feshbach resonances as 
demonstrated in recent scattering experiments [4, 13, 14]. 
Each resonance position marks an avoided crossing be¬ 
tween the atomic scattering threshold and a molecular 
bound state, which can be used to associate molecules 
from atom pairs [15]. 

We create and probe Er 2 dipolar molecules by us¬ 
ing standard magneto-association and imaging tech¬ 
niques [15]. Details of the production schemes are de¬ 
scribed in the Supplemental Material [16]. In brief, we 
begin with an ultracold sample of ^®®Er atoms in an op¬ 
tical dipole trap (ODT) in a crossed-beam configuration. 
The atoms are spin-polarized into the lowest Zeeman sub- 
level (j = 6, rrij = —6). Here, j is the atomic elec¬ 
tronic angular momentum quantum number and mj is 
its projection on the quantization axis along the magnetic 
field. To associate Er 2 molecules, we ramp the magnetic 
field across one of the low-field Feshbach resonances ob¬ 
served in Er [4, 13]. We experimentally optimize the 
ramping parameters, such as the ramp speed and the 
magnetic-field sweep interval, by maximizing the conver- 
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Figure 1. (color online). Er 2 weakly bound molecules, (a) 
Atom-loss spectrum [4] from 0 to 3 G and (b) near-threshold 
binding energy of the corresponding molecular states. The 
solid lines are hts to the experimental data and extrapolated 
to larger Ef, up to h x 500 kHz. The error bars are smaller 
than the symbols. 


sion efficiency. In our experiment we typically achieve a 
conversion efficiency of 15 %, which is a common value 
for boson-composed Feshbach molecules [15]. To obtain 
a pure molecular sample, we remove all the remaining 
atoms from the ODT by applying a resonant laser pulse. 
Our final molecular sample contains about 2 x 10^ Er 2 
Feshbach molecules at a temperature of 300 nK and at a 
density of about 8 x 10^^ cm“^ [16]. 

A central question regards the magnitude of the dipole 
moment owned by the magnetic molecules. We experi¬ 
mentally determine /r by using magnetic-field modulation 
spectroscopy, a technique which was successfully applied 
to alkali atoms [17-19]. With this method, we measure 
the molecular binding energy Ei, near the atomic thresh¬ 
old as a function of the magnetic field B. The bind¬ 
ing energy is related to the differential magnetic mo¬ 
ment of the molecules with respect to the atom-pair 
magnetic moment 2a. Here, /ia = —grrij^B = 6.98 
in the case of Er, where g = 1.16 is the Er atomic 
Lande factor. We thus extract g. by using the relation 
g = 2ga — \dEi,{B)/dB\. Our spectroscopic measurement 
begins with an ultracold atomic sample near a Feshbach 
resonance. We then add a small sinusoidal modulation 
to the bias magnetic field for 400 ms. The modulation 
frequency is varied at each experimental run. When it 
matches E^/h, prominent atom losses appear because of 
molecule formation. We trace the near-threshold molec¬ 
ular spectrum by repeating the measurement for various 
magnetic-field values. Figure 1 shows the Er 2 molecular 
spectrum in a magnetic field range up to 3 G. In our 
range of investigation, we identify four molecular energy 
levels, which near treshold exhibit a linear dependence on 
B. For each state, we obtain a different g value, ranging 
from 8 to 12 gs [20], as listed in Table I. 

For alkali-metal atoms, which possess much simpler 
interaction properties than lanthanides, theoretical ap¬ 
proaches based on coupled-channel calculations have 
been extremely successful in assigning the quantum num- 


Table I. Experimental and theoretical magnetic moments of 
four molecular states near the atomic threshold, Feshbach- 
resonance positions Bpr, dipolar lengths, outer turning points 
R* , and dominant quantum numbers i, J, and M. For conve¬ 
nience, the molecular states are labeled as gi with i = 1,... 4. 
The specified uncertainties correspond to the Icr statistical 
errors. 



Bfr 

m/m 

B 

ad 

R* 

\£,J,M) 



(G) 

Exp. 

Theo. 

(ao) 

(ao) 




0.91 

11.30(7) 

11.20 

1041(13) 

72.0 

|4,12, - 

12/-10/-9) 

M2 

2.16 

11.51(4) 

11.46 

1080(8) 

71.0 

|4,10, - 

10) 

M3 

2.44 

11.84(2) 

11.75 

1143(4) 

86.0 

|2,12, - 

10) 

M4 

2.47 

7.96(3) 

7.92 

517(4) 

57.0 

|6,10, - 

7/-6) 


bers of the molecular energy levels and reproducing 
molecular spectra [15]. However, a straightforward ex¬ 
tension of these methods to the lanthanide case is out 
of reach because of their complex scattering physics in¬ 
volving highly anisotropic interactions and many par¬ 
tial waves [13]. Inspired by work on alkali-metal colli¬ 
sions [21-24], we develop a new theoretical approach to 
identify the molecular quantum numbers, based on ap¬ 
proximate adiabatic potentials and on the experimentally 
measured g as input parameters. Our scattering model 
is detailed in the Supplemental Material [16], whereas we 
here summarize the central ideas of our approach. 

We first solve the eigenvalue problem of the full atom- 
atom interaction potential operator [16], whose eigenval¬ 
ues are the adiabatic potentials Un{R;B). The corre¬ 
sponding eigenfunctions read as \n;R) = J2i'^n,i{R)\i), 
where n = 1,2,... and Cn,i{R) are i?-dependent coeffi¬ 
cients. The molecular state jf) is uniquely determined by 
the set of angular momentum quantum numbers (£, J, 
M), where ^ is the molecular orbital quantum number, 
J = ji + 32 the total atomic angular momentum, and M 
its projection on the internuclear axis. 

To derive the corresponding “adiabatic” molecular 
magnetic moments, we calculate ~ —dUn{R] B)/dB 
at the position of the outer classical turning point R = 
R* . This choice is justified by the fact that most of the 
vibrational wavefunction is localized around R*. 

From the Hellmann-Feynman theorem it then follows 
that -^*1^(7?*)]^. Finally, we assume 

that for each Feshbach resonance a vibrational state is 
on resonance and we find the adiabatic potential that has 
a magnetic moment closest to the measured one within 
1 %. Once the best match is identified, the correspond¬ 
ing \n;R) sets the molecular state Ji), characterized by 
£, J, and M, with the largest, dominant contribution. 
In our range of investigation we observe d-, g-, and i- 
wave molecular states; see Table I. These states show 
several dominant M contributions. This fact is unusual 
and reflects the dominant role of the DDI, which cou¬ 
ples several adiabatic potentials and M components. As 
shown in Fig. 2, this mixing effect is particularly dom- 



















3 



Figure 2. (color online). Adiabatic magnetic moment as a 
function of magnetic field strength evaluated at the entrance 
channel energy. Each curve corresponds to the adiabatic mag¬ 
netic moment of one adiabatic potential Un{R; B). The mag¬ 
netic moments in the asymptotic limit of large B are given. 
The dashed vertical lines correspond to the field strength 
where we have observed Feshbach resonances. The red-filled 
circles represent the experimentally measured magnetic mo¬ 
ments at these resonance locations. 


inant below 10 G, where the DDI at R* is larger than 
the Zeeman interaction. Above 10 G, we predict /i to be 
equal to integer multiples of g^B [16]. 

As summarized in Table I, we find very good over¬ 
all agreement between the measured and the calculated 
molecular magnetic moments. For the largest observed /i, 
we calculate a corresponding dipolar length ad « 1150 ao. 
This value exceeds the typical range of the vdW poten¬ 
tials, setting the DDI as the dominant interaction in the 
system. Remarkably, Ud for Er 2 is comparable to the one 
realized with ground-state KRb molecules [25] , which are 
an extensively investigated case serving as a benchmark 
dipolar system. 

Following the methods introduced for KRb [26-28] , we 
test the dipolar character of Er 2 by performing scatter¬ 
ing experiments in a three- (3D) and in a quasi two- 
dimensional (Q2D) optical dipole trap. We control the 
DDI between molecules by tuning the dipole orientation, 
which is controlled by changing the direction of the mag¬ 
netic field and is represented by the angle 9 between the 
magnetic field axis and the gravity axis. Our experiment 
begins with the atomic sample trapped either in a 3D 
or in a Q2D ODT. The Q2D trap is created by super¬ 
imposing a vertically oriented, one-dimensional optical 
lattice [16]. After the magneto-association and the re¬ 
moval of the remaining atoms, we probe the number of 
molecules as a function of the holding time in the ODT. 
We perform measurements for the molecular states fii, 
/i 2 , and /i 4 [29]. For each of these states, we measure 
the collisional stability of the sample for both in-plane 
{6 = 90°) and out-of-plane {6 = 0°) dipole orienta- 



holding time (ms) 


Figures, (color online). Typical time evolution of the number 
of molecules for 8 — 90° (squares) and 9 — 0° (circles) in a 3D 

(a) and in a Q2D trap (b). The data refer to molecules in the 
state pi for the 3D case (a) and molecules in the state p 2 in 
Q2D (b). The insets in (b) show an illustration of molecules in 
pancake-shaped traps with out-of-plane (right) and in-plane 
(left) orientations. The solid lines are two-body decay fits to 
the data. The error bars for (a) and (b) are smaller than the 
data points and are not shown. The data points in (a) are 
obtained by averaging 5 independent measurements and in 

(b) about 50 measurements have been averaged. 

tion, and extract the corresponding relaxation rate co¬ 
efficients, j3± and /3||, using a standard two-body rate 
equation [30]. 

Figure 3 shows typical molecular decay curves in (a) 
3D and in (b) Q2D. In 3D, we confirm that the in¬ 
elastic decay does not depend on 9. We obtain /Jsd = 
1.3(2) X 10“^°cm^/s. This is a typical value for boson- 
composed Feshbach molecules, which undergo a rapid vi¬ 
brational quenching into lower-lying molecular states, as 
demonstrated with alkali atoms ]30]. Contrary, in Q2D 
the decay rates clearly depend on the dipole orientation. 
For each investigated molecular state, /3 _l is larger than 
/?ll. We find a reduction of losses of up to 30 % for out- 
of-plane orientation, for which the DDI is predominantly 
repulsive. The ratio increases with increasing p; 

see Table IF We note that stronger suppression of losses 
can be obtained using a tighter two-dimensional confine¬ 
ment [26], which is presently not reachable with our ex¬ 
perimental parameters. 

The reduction of losses in Q2D draws a natural anal¬ 
ogy with the observations obtained with KRb molecules 
[28]. From a comparative analysis between Er 2 and KRb, 
one can unveil universal behaviour attributed to the DDI, 
for systems being different in nature, but sharing a simi¬ 
lar degree of dipolarity. We thus theoretically study the 
scattering behavior of Er 2 using a theoretical approach 

















4 


Table II. Experimental and theoretical loss rate coefficients 
P for T = 400 nK and for various fi and 9 aX B = 200 mG. 
Uncertainties of /3 are statistical from fitting and systematic 
due to number density uncertainty. For the slightly different 
values of /i compared to Table I and the error discussion see 
Supplemental Material [16]. 



9^/lJ-B 

P± (10"® cm 
Exp. 

Vs) 

Theo. 

/3|| (10 ® cm 
Exp. 

Vs) 

Theo. 

9-i 

8.7(6) 

12.5±0.3±3.3 

6.00 

10.6±0.3±2.8 

4.79 

di 

10.9(5) 

9.5 ±0.2 ±2.5 

6.81 

7.3 ±0.1 ±2.1 

5.07 

d2 

11.7(3) 

11.3±0.2±2.9 

7.12 

8.6 ±0.2 ±2.3 

5.13 


similar to the one successfully applied to KRb. Our for¬ 
malism, which accounts for the DDI and the isotropic 
vdW interaction, is described in Refs. [16, 31]. 

We compute the Er 2 -I- Er 2 loss rate coefficients /3(T) 
in 3D and in Q2D for given values of /r, 0, and T. By 
averaging over a 3D and a 2D Maxwell-Boltzmann dis¬ 
tribution, we obtain the thermalized loss rate coefficients 
/3(T) in 3D and in Q2D, respectively. In 3D, we find 
a rate coefficient of 1.0 x 10“^°cm^/s at T = 300nK, 
which is close to the experimental value [32]. In Q2D, 
our calculations show that the collision dynamics at long 
range, and thus the value of /3, depends on the dipole 
orientation and monotonically increases with fi. As in 
the experiments, our calculations show that collisions for 
in-plane orientation {jPiP) lead to larger molecular losses 
than for out-of-plane orientation (/3||). In Table II, we 
compare theory and experiment. The absolute values of 
(3 agree within a factor of two. This difference is well 
explained by the fact that our model does not include 
details of the short-range physics, being the Er 4 poten¬ 
tial energy surfaces currently unknown [16]. 

Remarkably, the experimental and calculated ratios 
agree very well with each other; see Fig. 4. This 

suggests that for Er 2 Feshbach molecules is deter¬ 

mined by the DDI and not by the short-range physics, 
and that it can be correctly described using a point- 
like-dipole formalism [16]. Figure 4 shows the compara¬ 
tive analysis between bosonic '^^K®^Rb and ^®®Er 2 , and 
fermionic ^°K®^Rb and ^®^Er^®®Er based on our numer¬ 
ical calculations. Independent of the nature of the mag¬ 
netic or electric dipolar system, we find universal curves 
as a function of ad/a: one for bosons with a = Oho and 
one for fermions when d = Ovdw- Here, Uho is the har¬ 
monic oscillator length and Ovdw = is the 

vdW length with Ce the vdW coefficient. The faster in¬ 
crease of ^ for fermions with respect to bosons is due 
to the statistical fermionic suppression of ,5|| in Q2D that 
does not occur for bosons as explained in Ref. [33] . 

The universal behavior of ultracold dipolar scatter¬ 
ing has been previously pointed out in Ref. [34]. In the 
Wigner regime, we derive simple universal scaling laws 
for dipolar bosonic and fermionic molecules [16, 34]. For 



Figure 4. (color online). Universal loss rate ratio P±/P\\ as 
a function of ad/d for Er 2 (circles) and KRb (squares) for a 
fixed value Odn/aho = 4.85 corresponding to T = 400(40) nK 
and = 31.2(1) kHz [16]. The gray shaded area is due to 
the uncertainty of T. Here, d = Uho for bosonic molecules 
(filled symbols) and d = Ovdw for fermionic molecules (open 
symbols). The calculated loss rate ratios of Er 2 are compared 
with the experimental data for states ^i, /r 2 , and ^4 (trian¬ 
gles). 


bosons with ad, Oho > a^dw, which is the case of our Er 2 

molecules, we find ~ exp . 

> P||(r) V Oho / Oho Vaho/ 

For fermions with Od,avdw < Oho, ^ ( a w ) ' 

OdB = h/y/2T:mk^T is the thermal de Broglie wave¬ 
length. 


To conclude, our work reports on the study of strongly 
dipolar molecules created by pairing ultracold atoms with 
large magnetic dipole moments. We anticipate that 
our scheme can be generalized to other magnetic lan¬ 
thanide species and has the potential to open regimes 
of investigations, which have been unaccessible so far. 
First, the extraordinarily dense and rich molecular energy 
spectrum of Er opens the exciting prospect of cruising 
through molecular states of different magnetic moments 
or even creating molecular-state mixtures with dipole im¬ 
balance [19, 35, 36]. Second, in contrast to electric polar 
molecules where the electric dipole moment is zero in the 
absence of a polarizing electric field, magnetic dipolar 
molecules have a permanent dipole moment allowing to 
study the physics of unpolarized dipoles. In addition, 
strongly magnetic Feshbach molecules offer a novel case 
of study for scattering physics. These molecules are in 
fact diffuse in space with a typical size on the order of 
the vdW length. This novel situation can also have in¬ 
teresting consequences and trigger the development of 
extended scattering models, which account for multi¬ 
polar effects and truly four-body contributions when the 
molecule size becomes comparable to ad [37]. Finally, a 
very promising development will be to create fermionic 
Er 2 dipolar molecules where vibrational quenching pro¬ 
cesses are intrinsically suppressed because of the Pauli 
exclusion principle [38, 39]. 
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SUPPLEMENTAL MATERIAL 
Creation of Er 2 in 3D and Q2D 

We create Feshbach molecules using standard tech¬ 
niques of magneto-association across a Feshbach reso¬ 
nance. As demonstrated in Refs. [4, 13], Er features an 
enormous number of Feshbach resonances. Here, we fo¬ 
cus on the resonances observed below 3 G. In particu¬ 
lar, we first create an ultracold atomic sample of about 
3 X 10® ^®®Er atoms at a temperature of T « 150 uK, 
which is just above the onset of Bose condensation, see 
Ref. [4] . The atoms are confined into a three-dimensional 
(3D) crossed optical dipole trap with frequencies Vx = 
51.5(2) Hz, Vy = 13.2(3) Hz, and = 207(1) Hz. We 
choose magnetic fields of 1.4 G, 2.3 G, and 2.8 G for the 
molecular states /ri, ^ 2 , and ^ 4 , respectively. We then 
magneto-associate molecules by ramping the magnetic 
field 150 mG below the Feshbach resonance. The typi¬ 
cal ramp speed is 90mG/ms. After the molecule asso¬ 
ciation, we remove all the residual atoms from the op¬ 
tical dipole trap by applying a short laser pulse. The 
pulse is on resonance with the strong atomic transition 
at 401 nm [40] and has a duration of 1 fis with an intensity 
of ~ 40 mW/cm^. 

To realize a Q2D geometry, we superimpose a one di¬ 
mensional optical lattice beam to the system after finish¬ 
ing evaporation in the 3D trap. The lattice is realized 
from a retro-reflected laser beam at 1064 nm, propagat¬ 
ing along the vertical direction. The beam has a waist of 
250 /rm and a typical power of 8 W. As a result, the par¬ 
ticles are confined into an array of Q2D pancakes with 
frequencies = 33.0(3) Hz in the radial direction and 
= 31.2(1) kHz in the tightly confining axial direc¬ 
tion. We first load the lattice from the atomic sample 
and we then magneto-associate Er 2 in the lattice. The 
molecule conversion efficiency in the Q2D geometry is 
<5%, which is below the one observed in the 3D trap. 
With this scheme, we produce about 1.1 x 10^ molecules 
at a temperature of 400 nK, corresponding to a density of 
3.8 X 10^ cm“^. The molecules fill about 35 lattice layers. 

We control the molecular dipole orientation by chang¬ 
ing the orientation of the magnetic field. The orienta¬ 
tion is quantified in term of the angle 6 , which defines 
the angle between the quantization axis, set by the mag¬ 
netic held orientation, and the z-axis of the lattice trap. 
We prepare the molecular samples at either 0 = 0° or 
90°, correspondingly side-by-side (repulsive) or head-to- 
tail (attractive) dipolar collisions. The magnetic held is 
rotated by using three pairs of independently-controlled 
magnetic-held coils. We pay particular attention that 
when changing the orientation of the magnetic held we 
keep its magnitude constant. We check this by perform¬ 
ing radio-frequency spectroscopy between Zeeman sub- 
levels for different angles of rotation. We typically rotate 


the held within ^ 6 ms. 

For all our loss-rate measurements, we jump to a mag¬ 
netic held of about 200 mG after molecule association. At 
this held, Ej, is of the order of few hxl MHz. We choose 
to perform our measurement at this magnetic-held value 
because around 200 mG there are no Feshbach resonances 
and the molecular spectrum might be less dense. Using 
a Stern-Gerlach technique [19], we measure fi at 200 mG 
for all the three target molecular states. We hnd a slight 
shift of ^ in comparison with the values from Table H of 
a few percent to 10.9(5) i^b for /ii, 11.7(3) for /i 2 , and 
8.7(6) for ^ 4 . 

The given uncertainties for the measured loss rates in 
Table H are composed of a statistical error with one stan¬ 
dard deviation derived from htting a two-body rate equa¬ 
tion to the measured data, and a systematic uncertainty 
coming from number density calibration. Due to the dis¬ 
tribution of molecules across many lattice layers this is by 
far the greatest uncertainty in the Q2D geometry. The 
average 2D density and its uncertainty was calculated 
using a number-weighted average over occupied lattice 
layers similar to Ref. [28] . When calculating the loss rate 
ratio /3_l//3||, the systematic uncertainty in the density 
can be neglected as it is highly correlated for the mea¬ 
surement of j3± and /3||. 

Collision Formalism 

We briefly describe the theoretical formalism used in 
this article to determine the collisional properties of Er 2 
molecules in free space (3D collisions) and in an one¬ 
dimensional optical lattice (Q2D collisions), in an arbi¬ 
trary magnetic field B. More details can be found in 
Ref. [31, 33]. 

We use a time-independent quantum formalism based 
on spherical coordinates r = (r, Or^ipr) describing the rel¬ 
ative motion of two Er 2 molecules. The quantization axis 
z is chosen to be the confinement axis of the optical lat¬ 
tice. A spherical harmonic basis set, summed over dif¬ 
ferent partial waves ^ with projections mi on the quan¬ 
tization axis, is used to expand the total colliding wave 
function. The one dimensional optical lattice is supposed 
to be deep enough to consider the collision taking place 
in an individual pancake. One pancake is represented as 
an harmonic trap for the relative motion of reduced mass 

Wo — ^ (^) 

with uj = 27r^ and v = 31.2 kHz. The 3D collisions are 
recovered by setting 12 = 0. We consider molecules in the 
ground state of the harmonic oscillator. A given state of 
an Er 2 Feshbach molecule is described by a rather com¬ 
plicated linear combination of atomic states which cannot 
be precisely calculated as mentioned in the next section 
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of this Supplemental Material. Therefore we consider 
that the molecule has a magnetic moment of magnitude 
/i aligned along the magnetic field which makes an angle 
9 with the confinement axis. The interaction between 
two molecules is provided by the magnetic dipole-dipole 
interaction 


Kid 


(1 — 3cos^(0r. — 9)) 


( 2 ) 


We also used an isotropic Er 2 -I- Er 2 van der Waals in¬ 
teraction given by 



with Ce = 4 X 1760 = 7040 a.u. which amounts to four 
times the value of an isotropic atom-atom coefficient of 
1760 a.u. from the theoretical work of Ref. [12] . Note that 
an alternative value of 1723 a.u. based on observed transi¬ 
tions was obtained in Ref. [13]. The Schrddinger equation 
is solved for each radial intermolecular separations r us¬ 
ing a log-derivative propagation method. Matching the 
colliding wavefunction and its derivative with appropri¬ 
ate two-dimensional asymptotic boundary conditions at 
long-range [33] provides the cross section and the rate co¬ 
efficient as a function of the collision energy for any arbi¬ 
trary configurations of magnetic fields and confinements. 
Averaging the cross sections over a 3D and 2D dimen¬ 
sional Maxwell-Boltzmann distribution provides the cor¬ 
responding thermalized rate coefficients /3(T) for a given 
temperature. 

At short range, we assume that the molecules undergo 
a full loss mechanism process with a unit probability (it 
can be either an inelastic or a possible reactive process). 
This assumption, which corresponds to the so-called uni¬ 
versal regime in ultracold collisions, considers that the 
physics is independent of the initial short-range scatter¬ 
ing phase-shift [41] of the full potential energy surfaces 
of Er 4 . This is what it is usually assumed for theory as 
nothing is known about this potential energy surface at 
short range. Then, if the magnitude of the rates differs 
between experiment and theory, one can learn that an 
experimental system deviates from this universal regime 
and short-range effects play a role. 

To circumvent this, it is more convenient to compute 
the ratio of the theoretical rates of two different mag¬ 
netic field orientations since we will start with the same 
short-range physics condition for both orientations, and 
compare it with the corresponding experimental ratio. 
An analysis based on the universal behavior of dipolar 
collisions in confinement of Ref. [34] using a Quantum 
Threshold model leads to the following formula for the 
ratio /3j_(r)/,d||(T). Eor bosons, using Eq. 30 of Ref. [34] 
to describe j3±_ (dipole dominated) and Eq. 32 of the same 
reference for /3|| (confinement dominated) we find 


PAT) 

P\\iT) 



^ho 


(4) 


when Od, Uho > Ovdw for a fixed value of ade/aho where 
OdB is the thermal de Broglie wavelength. Eor fermions, 
using Eq. 16 of Ref. [34] to describe /3_l (dipole domi¬ 
nated) and Eq. 14 of the same reference for /Sy (van der 
Waals dominated), along with Eq. 27, we find 


A 

P\\ 




fer 

V flvdW J 


(5) 


when ad, Ovdw < Oho- These formulas suggest to plot the 
ratio as a function of Od/oho for bosons for a fixed ratio 
OdB/oho = v/ksT and as a function of Od/ovdw for 
fermions, as it has been done in Fig. 4 for the magnetic 
dipolar molecules of Er 2 and the electric polar molecules 
of KRb. 


Adiabatic Model 

In Ref. [13] we presented the theoretical bosonic- 
erbium Feshbach spectra derived from coupled-channels 
calculations. We concluded there that such first-principle 
evaluations can not quantitatively capture the complex 
scattering behavior of Er. In fact with the current com¬ 
puting capabilities, the calculations can not be converged 
with respect to the number of basis states required to ex¬ 
plain the experimental Feshbach-resonance density. For 
this reason, we developed a novel approach based on adi¬ 
abatic potentials (adiabats) Un{R;B). 

Our adiabatic model starts from the Hamiltonian H = 
— {tP/2mr)(P/dE? + K(i?). The first term is the radial 
kinetic-energy operator with R describing the orienta¬ 
tion and the separation between the two atomic dipoles, 
and nir is the reduced mass. The second term of the 
Hamiltonian is the potential operator V{R), which de¬ 
scribes the Zeeman and interatomic interactions. It reads 
V{R) = /{2mrR'^) + Hz + W^^^^{R) and incorporates 

the rotational enerj;y operator with molecular orbital an¬ 
gular momentum £, the Zeeman interaction of two atoms 
Hz, and the electronic potential operator be¬ 

tween the particles. Our model assumes that the relative 
vibrational motion of two Er atoms is slow compared to 
the timescales of the rotational, Zeeman, and “electronic” 
atom-atom interactions. 

The Zeeman interaction is Hz = giJ^B{,jiz + j 2 z)B. 
Here, g = 1.16 is the Er g-factor, a magnetic field B 
is aligned along the z direction, and jA is the z com¬ 
ponent of the angular momentum operator jl of atom 
i = 1,2. The electronic potential operator W®'®“(.R), de¬ 
scribed in Refs. [11, 13, 42], is anisotropic, as it depends 
on the orientation of R. At large separation R, lT®*®®(f?) 
is given by the magnetic dipole-dipole interaction plus 
both the isotropic and anisotropic contribution of the 
van der Waals interaction. For i? —>■ oo the interaction 

^el®c(^) ^ 

The Hamiltonian is evaluated in the basis jz) = 
\Uih)JM)Yg^^{R), where J = ji + j 2 and Y(mt{R) is 








a spherical harmonic. It conserves mi + M and parity 
p = (—1)^. In addition, for bosonic isotopes (—1)^+'^ = 1. 
We focus on ultracold collisions between atomic states 
IjiTOi) = \j2'm2) = |6, — 6) and, therefore, only include 
basis functions satisfying mi + M = —12. We limit the 
included partial waves to even ^ < 6 and thus to states 
with even J, as the “adiabatic” magnetic moments of the 
resonances quickly converge with the included number of 
partial waves (In our calculation there is one s-wave chan¬ 
nel, four d-wave channels, nine g-wave channels, and 16 
i-wave channels.). 

The adiabats Un{R; B) with n = 1,2,... are eigen¬ 
values of the operator V{R) at a given field strength 
B. Their eigenfunctions are |n;i?) = Cn^i{R)\i} with 
i?-dependent coefficients Cn,i{R)- Note that we neglect 
the coupling between Un{R] B) due to the radial part of 
kinetic-energy operator. 

Figure SI shows the adiabats at B = 2.44 G. The scat¬ 
tering starts from the s-wave entrance channel correlating 
to the energetically lowest adiabat. All other potentials 
either have a centrifugal barrier and dissociate to two 
atoms with M = —12, or dissociate to closed-channel 
Zeeman sublevels with M > —12. We distinguish four 
groups of potentials, each associated with a dominant 
partial wave £. Within a group, the potentials are split 
by the Zeeman energy and the magnetic DDI and disso¬ 
ciate at different atomic thresholds. For each potential 
Un{R]B) we can further assign the dominant J and M, 
where T = Ji + j 2 is the sum of electronic angular mo¬ 
menta of two atoms and M is the projection of J on 
the internuclear axis. The figure also shows an exam¬ 
ple of predominantly d-wave Feshbach molecules with an 
outer classical turning point R*. Its “adiabatic” molec¬ 
ular magnetic moment is to good approximation given 
by « —dUn{R*',B)/dB, where we further use that 
most of the vibrational wavefunction is localized around 
R*. Interestingly, we observe that the value quickly 
converges with the number of included I (even £ < 6 is 
sufficient) and that it strongly depends on the DDI but 
only weakly on the vdW dispersion potential. In fact, at 
R* the DDI dominates over the anisotropic part of the 
dispersion potential. 

The adiabatic magnetic moment of a resonance is 
given by = -dEnv{B)/dB « -dUn{R*\B)/dB, 

where we realize that to good approximation most 
of the adiabatic vibrational wavefunction is localized 
around the outer classical turning point. We further 
note that dUn{R*B)/dB = {n; R*\dHz/dB\n; R*) from 
the Hellmann-Feyman theorem and, hence, = 

—gpB AIi|ci(i?*)p, where Mi is the total atomic pro¬ 
jection quantum number of state |i). We assign a res¬ 
onance by the quantum numbers of the basis state |i) 
for which |ci(i?*)p is largest and note that the absolute 
value of the magnetic moment of a resonance is always 


smaller that VlgpE ~ 14/i_B- 
(a) (b) 



Figure SI. Adiabatic interaction potentials of two Er atoms 
at medium- (a) and long- (b) interatomic separation R. The 
calculation is performed at B = 2.44 G with mi -\- M = —12 
and includes only states with even I < &. The zero of en¬ 
ergy is at the dissociation limit of two \ j, m) = |6, —6) atoms. 
Black, green, red, and blue curves indicate the dominant £- 
wave character. The horizontal black line indicates a d-wave 
Feshbach molecule with an outer turning point R = R* reso¬ 
nant with the s-wave entrance channel. Panel (b) also shows 
the M projection for each of the Zeeman dissociation limit. 

We further assume that the non-adiabatic coupling 
between the adiabatic potentials is significantly smaller 
than their spacings for R < lOOoo. Then a weakly bound 
level of adiabatic potential n can lead to a Feshbach res¬ 
onance when its energy Env{B) coincides with the en¬ 
trance channel energy. The outer turning point R* of 
this level satisfies Un{R; B) = 0. The resonance acquires 
a width due to non-adiabatic coupling to the entrance 
channel. 

Finally, we determine the approximate quantum num¬ 
bers of experimentally-observed resonances with < 3 
G, listed in Table I, based on a comparison of the exper¬ 
imental magnetic moment with those predicted by the 
adiabatic model at the same resonant field. We find that 
for these resonances there exist adiabats with a magnetic 
moment that agrees within 1% uncertainty with the ex¬ 
perimental values. A study of the largest coefficients 
Cn^i{R) at R = R* then enables us to assign the dom¬ 
inant quantum states shown in Table I. 

Figure 2 (main text) shows the magnetic-field de¬ 
pendence of the adiabatic magnetic moment at the en¬ 
trance channel energy for each of the adiabatic poten¬ 
tials Un{R;B). We see that for i? > 10 G the mag¬ 
netic moment values equal integer multiples of ggs cor¬ 
responding to those of the atomic limits. For smaller field 
strengths the adiabatic magnetic moments show mixing 
of the Zeeman sublevels. Here, the magnetic moment 
value depends on the magnetic dipole-dipole interaction 
but only weakly on the strength and anisotropy of the 
dispersion potential. The figure also shows our experi¬ 
mentally studied Feshbach resonance locations as well as 
their magnetic moments see Table I. 



















